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1.  Objectives:  This  project  is  developing  mathematical  foundations  for  sensor 
management  systems.  Sensor  management  can  be  viewed  as  the  process  of  selecting 
sensors,  sensor  modes  cind  sensor  dwell  points  to  optimize  the  probability  density 
characterizing  a  set  of  targets  in  a  region.  Thus,  there  are  two  parts  to  this  problem.  In 
order  to  support  this  approach,  representations  of  multitarget  probability  densities  must 
be  developed  that  are  rich  enough  to  model  the  uncertainty  between  quantities  such  as  the 
number  of  targets  in  a  region,  their  locations  and  target  classes.  Then,  given  this 
probability  density,  measures  of  sensor  effectiveness  must  be  computed  and  optimized. 
This  work  is  a  collaboration  between  Lockheed-Martin  Tactical  Defense  Systems,  the 
Minnesota  Center  for  Industrial  Mathematics  (MCIM)  and  Mr.  S.  Musick  at  the  USAF 
Wright  Laboratory. 

2.  Status  of  Effort:  Sensor  management  requires  optimizing  sensor  effectiveness  across 
multiple  conflicting  objectives  such  as  a  desire  to  localize  and  classify  known  targets 
while  attempting  to  detect  new  ones  as  well.  This  requires  joint  representation  of  the 
probability  density  for  the  number,  location  and  class  of  a  collection  of  targets.  In  order 
to  solve  this  problem  and  to  study  it  in  a  simple  setting,  we  developed  the  notion  of  Joint 
Multitarget  Probabilities  for  detection  and  tracking.  During  the  past  year  this  was 
extended  to  the  problem  of  sensor  management  with  target  classification. 

A  related  issue  that  arises  in  this  context  is  determining  how  these  probability  densities 
evolve  in  time.  The  resulting  densities  are  often  very  complicated.  Even  in  single  target 
cases,  determining  the  time  evolution  of  these  densities  is  a  challenge  requiring  new 
algorithmic  approaches.  In  certain  cases  this  evolution  is  characterized  by  a  partial 
differential  equation  called  the  Fokker-Planck  equation  (FPE)  and  its  generalizations. 
Although  this  has  long  been  known,  solution  of  the  FPE  has  generally  been  regarded  as 
unfeasible  for  real-time  applications.  Now,  recent  advances  in  computers  and  numerical 
methods  may  have  changed  this.  To  explore  the  viability  of  real-time  FPE  solutions,  we 
developed  several  solution  schemes  and  applied  them  to  test  problems  with  very 
promising  results. 

Some  of  our  past  work  on  sensor  management  for  multitarget/multisensor  developed 
discrimination  gain  for  the  case  of  a  collection  of  targets  described  by  conventional 
Kalman  filters.  In  a  generalization  of  this  work,  a  method  for  computing  the  expected 
discrimination  gain  was  developed  for  multiple  model  Kalman  filters.  This  is  being 
applied  in  a  more  realistic  setting  under  a  transitioning  Small  Business  Innovative 
Research  program. 

3.  Accomplishments:  This  year's  accomplishments  are  in  three  areas:  Extension  of  a 
Kalman-filter  based  discrimination  metric  for  sensor  management  to  a  more  general  class 
of  filters;  extension  of  sensor  management  based  on  Joint  Multitarget  Probabilities  to 
incorporate  multiple  sensor  modes  and  target  classification;  and  the  development  of  fast 


methods  to  solve  the  Fokker-Planck  equation  for  real-time  non-linear  filtering 
applications. 

Expected  discrimination  gain  is  a  useful  general  measure  of  sensor  effectiveness  that  was 
developed  for  targets  being  tracked  by  Kalman  filters.  Real  applications  often  require 
more  sophisticated  filters  which  characterize  target  motion  in  terms  of  both  continuous 
dynamical  states  and  discrete  states  referred  to  as  dynamics  modes.  A  measure  of 
discrimination  gain  was  derived  to  apply  to  this  case.  In  this  case,  there  are  two  factors 
that  contribute  to  the  expected  discrimination  gain:  the  reduction  in  uncertainty  about  the 
target's  kinematics;  and  the  reduction  in  the  uncertainty  about  the  target's  flight  mode,  e.g. 
whether  the  target  is  turning  or  in  straight  and  level  flight.  A  conference  paper  has  been 
accepted  on  this  work  [7]  and  a  more  detailed  journal  paper  has  been  submitted  [2]. 

Earlier  work  on  discrimination  gain  for  static  detection  problems  [1]  was  generalized  to 
the  problem  of  dynamic  targets  and  dynamic  targets  with  classification  and  sensors  with 
multiple  modes  to  select  from  [8].  In  this  case,  discrimination  can  be  optimized  by 
determining  where  to  direct  the  sensors  and  which  sensor  mode  to  use.  For  example,  in  a 
test  problem  with  a  sensor  that  can  either  detect  or  classify  targets,  discrimination  is 
initially  optimized  by  using  the  sensor  in  its  detection  mode.  After  the  targets  have  been 
detected  and  localized  the  classification  mode  begins  to  have  high  discrimination  value. 
This  method  is  very  compute  intensive.  In  collaboration  with  Mr.  S.  Musick,  a 
simplifying  scheme  based  on  a  product  representation  of  the  joint  density  is  being 
developed  [9]. 

In  order  to  compute  the  expected  discrimination  gain  when  a  sensor  mode  and  dwell  point 
is  selected,  the  time  evolution  of  the  probability  density  between  measurements  must  be 
computed.  This  can  be  determined  by  the  Fokker-Planck  equation  FPE  (or  closely  related 
equations).  In  collaboration  with  researchers  at  the  U.  of  MN  we  explored  several 
methods  for  solving  the  FPE  in  real  time.  With  Dr.  M.  Kouritzin,  a  convolution-based 
scheme  for  altitude  tracking  was  developed  and  evaluated  [10].  Several  techniques  for 
solution  of  the  FPE  using  finite  difference  schemes  were  explored.  As  part  of  his 
Master's  degree  work  at  the  Minn.  Center  for  Industrial  Mathematics,  Mr.  A.  Zatezalo 
developed  and  tested  four  FPE  solvers  against  a  model  problem  typifying  tracking 
applications.  Dr.  Kastella  supervised  this  work  while  Mr.  Zatezalo  interned  with 
LMTDS-E  in  the  summer  of  '96.  Dr.  N.  Krylov  is  Mr.  Zatezalo's  advisor.  Valuable 
insight  into  the  proper  way  to  treat  the  boundary  condition  for  this  type  of  problem  was 
provided  by  Prof  Friedman.  The  four  FPE  solvers  we  studied  [11]  were:  the  Alternating 
Direction  Implicit  (ADI)  scheme;  the  method  of  fractional  steps;  an  Euler  explicit 
scheme;  and  a  multigrid  scheme.  For  this  problem  ADI  and  the  method  of  fractional 
steps  were  found  to  be  about  two  orders  of  magnitude  more  efficient  than  the  Euler 
explicit  and  multigrid  schemes. 

Based  on  this  study,  the  Alternating  Direction  Implicit  (ADI)  scheme  seems  very 
promising  for  sensor  management  and  non-linear  filtering  applications.  ADI  is  a  so- 
called  0(N)  algorithm,  which  means  that  its  computational  complexity  grows  linearly 


with  the  size  of  the  finite-element  grid  used  to  discretize  the  problem.  In  an  effort  to 
assess  the  feasibility  of  ADI-based  non-linear  filtering,  a  diverse  set  of  model 
applications  have  been  formulated.  These  model  applications  are:  1)  low-SNR  image 
tracking  with  nonlinear  target  dynamics  and  non-Gaussian  measurements;  2)  missile 
altitude  tracking  with  severe  multipath  interference  and;  3)  joint  tracking  and 
identification  using  RADAR  position  measurements  and  range-profile  measurements  (an 
application  suggested  by  Mr.  Musick).  Two  conference  papers  on  this  work  have  been 
accepted  for  publication  [11, 12],  one  journal  paper  has  been  submitted  [3]  and  several 
more  are  nearing  completion  [4-6].  In  all  of  these  applications,  NLF  significantly 
improvements  estimation  performance  at  a  supportable  cost  in  computational  load.  For 
example,  in  low-signal/noise  ratio  image  tracking,  it  provides  an  order  of  magnitude 
improvement  in  performance  against  maneuvering  targets.  In  simulation  tests  for 
tracking  in  multipath  interference,  NLF  successfully  tracks  targets  that  are  essentially  un- 
trackable  by  other  means.  In  the  joint  tracking  and  recognition  application,  NLF  reduces 
the  time  to  classify  the  target  by  about  50%,  compared  to  a  direct  maximum-likelihood 
classifier. 
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6.2  Transitions: 

An  algorithm  using  discrimination-based  sensor  management  with  Interacting  Multiple 
Model  Kalman  Filters  (IMMKF)  is  currently  being  incorporated  into  a  Phase  I  Small 
Business  Innovative  Research  (SBIR)  contract  called  "Sensor  Management  Across 
Multiple  Platforms  (SMAMPS)". 
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